home *** CD-ROM | disk | FTP | other *** search
/ Amiga Format CD 46 / Amiga Format CD46 (1999-10-20)(Future Publishing)(GB)[!][issue 1999-12].iso / -in_the_mag- / reader_requests / scilab / demos / ode / lorenz.dem next >
Text File  |  1999-09-16  |  347b  |  15 lines
  1. deff('[ydot]=lorenz(t,y)',...
  2. "x=y(1);...
  3. a=[-10,10,0;28,-1,-x;0,x,-8/3];...
  4. ydot=a*y")
  5. deff('[j]=jacobian(t,y)',...
  6. "x=y(1);yy=y(2);z=y(3);...
  7. j=[-10,10,0;28-z,-1,-x;-yy,x,-8/3]")
  8. comp(lorenz);comp(jacobian);
  9. y0=[-3;-6;12];t0=0;step=0.01;t1=10;
  10. instants=t0:step:t1;
  11. y=ode(y0,t0,instants,lorenz,jacobian);
  12. xbasc(0);param3d(y(1,:),y(2,:),y(3,:))
  13.  
  14.  
  15.